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form factors to the continuum limit and to the physical values of the light-quark masses. 
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1. Introduction 



Flavor-changing neutral-current b — > s decays play an important role in constraining models 
of new physics. In addition to the widely-studied mesonic decays B — >■ K*y and B — >■ K [ *h + r, 
baryonic b ^ s decays decays such as Aj, — > A/ and A/, — > Al + £~ are also worth investigating 
since the spin of the baryons allows the construction of observables that are sensitive to the helicity 
structure of the effective weak Hamiltonian [1]. The CDF collaboration has recently observed the 
decay A/, — > Ap. + pi~ for the first time [2], and new results are expected from LHCb. 

To calculate the decay amplitudes for A& — > Ay and A& — > A£ + £~ , the hadronic matrix elements 
(A(p',j')| sTb \Ab(p,s)) need to be determined for F G {7^, 7^75, q v o )lv , q v o^ v Ys} (where q = 
p — p'), resulting in ten independent form factors. The situation simplifies when using heavy-quark 
effective theory (HQET) for the b quark. At leading order in heavy-quark effective theory, only 
two independent form factors remain, and one has [3] 

{A{ V \s')\sTQ\A Q {v,s)) =u( V ' ,s')[F Y (p' ■v) + fF 2 (p' ■ v)] r«T(v,j). (1.1) 

Here, v is the four- velocity of the Aq, and the form factors F\ and F2 are functions of p' ■ v, which is 
equal to the energy of the A baryon in the Ag rest frame. In our analysis, it proves more convenient 
to work with the linear combinations 



F+=F 1 +F 2 , F-=Fi-F 2 , (1.2) 

instead of F\ and F%. In the following, we report on our calculation of these two form factors using 
lattice QCD. For the static heavy quark Q, we set v = (1,0,0,0), and we use a lattice HQET action 
with one iteration of HYP smearing for the gauge link in the time derivative [4]. For the up, down, 
and strange quarks, we use a domain- wall action [5]. Our calculations are based on the 2+1 flavor 
RBC/UKQCD gauge field ensembles described in Ref . [6] . 



2. Extracting the form factors from correlation functions 



In our two-point and three -point correlation functions, we use the following interpolating fields 
for the Aq and A baryons, 



AQ a = e abc (C^)p 7 d a p u r Q c a , 



'a ~b /-ic 



A a = e abc (CY 5 )p 7 u a p d° r s c a 



(2.1) 



where the tilde on the up, down, and strange-quark fields indicates gauge-covariant Gaussian smear- 
ing. In the three-point functions, we use an ^(a)-improved discretization of the continuum HQET 
current, which is given by [7] 



r (HQET) 



(mi,) = U(ntb,a l ) 2f 



QTs + cp a) 



i-W h ) 



QTs + c { Y Psa) aQTy-Vs 



(2.2) 



The matching factor 5° and the improvement coefficients Cp and Cp have been computed in 



one-loop lattice perturbation theory in Ref. [7]. The factor U (nib, a ) provides two-loop running 
in continuum HQET from = a -1 to fi = ntb- 
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Figure 1: Propagator contractions for cfl(V, p',f ,t') (left) and C^ bw) (r, p',f,f - f') (right). The vertical 
thick lines indicate the static heavy-quark propagators. 



We compute "forward" and "backward" three-point functions originating from a common 
source point (xo,x), 

cg>(T, p',/,0 = -£e-^'^-yUA s (x ,x) 4 BQET) \x Q -t + t',y) A Qa (x - t,y)) , (2.3) 



c% bw \r, P ',t,t-t') = X^ , ' p '- (y " x) ( A e«(^+^y)^ QET) (^o+? , ,y) A 5 (*o,x)). 



(2.4) 



As is apparent from Fig. 1 , the three-point functions do not require sequential domain-wall propa- 
gators and can be computed efficiently for arbitrary values of t, t', F, and p'. We then construct the 
ratio 

4 Tr [c® (r, p', t, t') d 3 < bw > (r,p',t,t- 1')] 



(2.5) 



Tr[C( 2 A)(p',f)]Tr[C( 2 ' A f)(0] 

where C^ 2 ' A ^ (p', t) and C^ 2,Aq ^ (f) are the A and Aq two-point functions. By inserting complete sets 
of states and using Eq. (1.1), one finds that, for T equal to any product of y^'s, 



( E A + m A 



-A 



[F + ] 2 + ..., if[r,7°]=o, 
[f_] 2 + ..., if {r,r°} = o. 



(2.6) 



Here, F± are the form factors defined in Eq. (1.2), and the ellipsis indicates excited-state contribu- 
tions that decay exponentially with the time separations. To increase statistics, we average the ratio 
over multiple gamma matrices and define 

a+tfjj) = \ [^(i, P ',t/)+^(ff^,t/)+£g(fr l ,p',t/)+^(r l f^,t/)] , (2.7) 



(2.8) 



(Replacing T by Tyo would not give new information because JqQ = Q.) For a given value of 
|p'| 2 , we then average M±(\)i ,t,t') over the direction of p', and we denote these direction-averaged 
quantities as &±(\p'\ 2 ,t,t'). Because of the symmetric form of our ratio (2.5), at a given source- 
sink separation t, the excited-state contamination will be smallest at t' = t/2. We therefore define 
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Set 


A? 




(sea) 

am s 


sea) 

am u 4 


a (fm) 


(val) 
am s 


(val) 
am u,d 


(vv) 

mk 


(vv) 


C14 


24 3 


x64 


0.04 


0.005 


0.1119(17) 


0.04 


0.001 


245(4) 


761(12) 


C24 


24 3 


x64 


0.04 


0.005 


0.1119(17) 


0.04 


0.002 


270(4) 


761(12) 


C54 


24 3 


x64 


0.04 


0.005 


0.1119(17) 


0.04 


0.005 


336(5) 


761(12) 


C53 


24 3 


x64 


0.04 


0.005 


0.1119(17) 


0.03 


0.005 


336(5) 


665(10) 


F23 


32 3 


x64 


0.03 


0.004 


0.0849(12) 


0.03 


0.002 


227(3) 


747(10) 


F43 


32 3 


x64 


0.03 


0.004 


0.0849(12) 


0.03 


0.004 


295(4) 


747(10) 


F63 


32 3 


x64 


0.03 


0.006 


0.0848(17) 


0.03 


0.006 


352(7) 


749(14) 



Table 1: Properties of the gauge-field ensembles and propagators. Here, m n and m^ (both in units 
of MeV) are the masses of the pion and the pseudoscalar ss meson (without disconnected contributions), 
corresponding to the valence quark masses anv^' and am{ val \ respectively. 

the new quantities 

R + (\p'\ 2 ,t) = J-^^ + (\p'\ 2 ,t,t/2), (2.9) 

V t A + m A 

MIp'IV) = J ir ^0-(\V\ 2 ,t,t/2), (2.10) 

V E A - m A 

which, according to Eq. (2.6), are equal to the form factors F±(v • p') up to excited-state effects that 
decay exponentially with the source-sink separation t. We obtain the energy FaGp'I 2 ) and mass 
m A appearing in Eqs. (2.9) and (2. 10) from fits to the A two-point functions for the same data set. 

3. Data analysis 

Our analysis uses seven different data sets C14, F63 with parameters as shown in Table 1. 
We computed the ratios (2.7) and (2.8) for all possible lattice momenta up to |p'| 2 = 9 • (2n) 2 /L 
and for all source-sink separations in the range from 4 <t /a < 15 at the coarse lattice spacing and 
5 < t/a < 20 at the fine lattice spacing. Sample results for &±(\p'\ 2 ,t,t f ) are shown in Fig. 2 as 
a function of the current-insertion time t'. Note that there are plateaus in t', but the non-negligible 
dependence on t indicates that there are still excited state-contributions. This can be seen more 
clearly in Fig. 3, where the corresponding results for the quantities (2.9) and (2.10) are plotted 
against t. To isolate the ground-state contributions (i.e., the form factors F±), we extrapolate R± to 
t = oo, allowing for excited states using the ansatz 

Bg(t )=F j* +A g exp[ _ S £ t]f ( 3.i) 

where the label i = C14, ...,F63 denotes the data set, and n labels the momentum squared: |p'| 2 = 
n • (2n) 2 /L 2 . Because the energy gap parameters d l ± in Eq. (3.1) are positive by definition, we 
rewrite them as 8± = exp(/±") • (1 GeV), and use l l £ along with F±" and A±" as the fit parameters. 
Where necessary, we exclude a few points at the smallest t from the fits to get good # 2 /d.o.f. At 
fixed momentum n, we perform the fits simultaneously for the different data sets i. Since 2n/L 
(in GeV) is equal within uncertainties for the coarse and fine lattice spacings, we know (from prior 
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Figure 2: Example results for the ratios i%±(\p\ ,t,r) from the C54 data set. The data shown here are for 
|p'| 2 = 4 • (2n/L) 2 , at source-sink separations (from left to right) t/a = 6,8, 10. 




1.6 




Figure 3: Example results for R±(\p'\ 2 ,t), along with fits using Eq. (3.1). Left panel: C54 data set; right 
panel: F43 data set. 



studies of the hadron spectrum on the lattice) that the physical energy gaps must be of similar size 
for the different data sets i. With this knowledge, we use Bayesian constraints that limit differences 
|Z±" — l± n \ to reasonable values (details will be given in Ref. [8]), which improves the stability of 
the fits. 

Having performed these fits, we noted that at each momentum n and data set i, the energy gap 
parameters l'_f and returned from the fit were equal within uncertainties. This is expected as 
long as the relevant excited states have non-zero matrix elements in both R + and R . We therefore 
performed new fits with shared energy gap parameters if = if = l l,n . These new, coupled fits 
had values of % /d.o.f as good as or better than the separate fits, and the extracted form factors 
were consistent with those from separate fits (see Fig. 4), so we use F±" from the coupled fits in 
the further analysis. To estimate the systematic uncertainty resulting from the choice of the f m j n 's, 
we compute the shifts in F±'" when increasing all ? m i n 's by one unit, and we add these shifts in 
quadrature to the statistical uncertainties. 

The last step of the analysis is to interpolate the dependence of the form factors on p' ■ v = E\ 
using a suitable smooth function, and extrapolate to the continuum limit and the physical u, d, and 
s-quark masses. Since E\ depends strongly on the quark masses, it is better to consider the form 
factors on the lattice as functions of E\ — M^. At the present level of statistical precision, and for 
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Figure 4: Fit results for F_£" from the data sets i = C54 (left) and i = F43 (right). Open symbols show results 
obtained using Eq. (3.1) with separate energy gap parameters If and if. Filled symbols show the results 
from coupled fits with If = If = l' n . The points are offset horizontally for clarity. 



Parameter Result 

N+ 3.188±0.268GeV 2 

X + 1 .852 ± 0.074 GeV 

N- 4. 124 ±0.750 GeV 2 

X_ 1.634 ±0.144 GeV 



Table 2: Preliminary results for N± and X± from fits using Eq. (3.2). 



the energy range considered here, we find that generalized dipole fits of the form 

N± 



i.n 

A 



in 



i \2 



X±+c L ± 



K) 2 - K J 



Phys\2 



+ c s + ■ 



K/-(< ys ) 2 



(3.2) 
(3.3) 



with parameters N±, X±, d±, c/ ±, and c Vj ± describe our data well (# 2 /dof = 0.59 for F + , and 
% 2 /dof = 0.43 for F_). In the continuum limit and for the physical light and strange-quark masses 
(we use m^ ys = 686 MeV [9]), these functions reduce to F± = N±/(X± +E\ — m\) 2 , which contain 
only the fit parameters N± and X±. Our preliminary results for these parameters are given in Table 
2. Plots of the fitted functions are shown in Fig. 5. 



4. Outlook 



We have performed the first lattice calculation of the Aq — > A form factors F + = F\ + F2 and 
F = F\ — F%. Using a ratio technique with a wide range of source-sink separations, we have 
achieved a high level of statistical precision. The dominant systematic uncertainties in our results 
are associated with the use of one-loop perturbative current matching, finite-volume effects, the 
naive linear extrapolations in the light-quark masses, and the continuum extrapolations. We es- 
timate that the total systematic uncertainty is below 10%; more details will be given in Ref. [8]. 
There, we will also present results for the differential branching fraction of — > ApL + pL~ . 
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Figure 5: Preliminary fits of the form factor data for F + and F using Eq. (3.2). The fits includes all seven 
data sets (see Table 1), but only three data sets are shown to save space. The bottom-right plot shows the 
fitted functions evaluated in the continuum limit and at the physical values of the light and strange-quark 
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